Entanglement of electron spins in superconductors 
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We investigate entanglement of two electron spins forming Cooper pairs in an s-wave supercon- 
ductor. The two-electron space-spin density matrix is obtained from the BCS ground state using a 
two-particle Green's function. It is demonstrated that a two spin state is not given by a spin singlet 
state but by a Werner state. It is found that the entanglement length, within which two spins are 
entangled, is not the order of the coherence length but the order of the Fermi wave length. 
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Entanglement, referring to the nonlocal quantum cor- 
relation between subsystems, is one of the key resources 
in quantum telcportation, quantum communication, and 
quantum computation JJ. A study of entanglement in 
many-body systems is of importance for not only its ap- 
plication to quantum information processing but also giv- 
ing us new insights on its relevant physics. For example, 
entanglement in many spin systems has been investigated 
in connection with quantum phase transition [2|. For 
a non-interacting electron gas, the entanglement length 
within which two electron spins are entangled is the order 
of the Fermi wave length Xp = 2iT/kF 0, 3 . 

A solid state entangler, analogous to the parametric 
down conversion producing a pair of entangled photons 
in quantum optics, is a device to generate a pair of en- 
tangled electrons in controlled way. It is of interest for 
the realization of scalable solid-state quantum comput- 
ers. There have been various proposals to create entan- 
gled pairs in solid state systems; spin-entanglement via a 
quantum dot .Si, lq| or via a magnetic impurity pi , gener- 
ation of entangled electron spins by extracting a Cooper 
pair out of a superconductor H i, [lOl [H [H El 111 , 
entangled electron-hole pairs in a degenerate electron 
gas [ij, [13 , two particle orbital entanglement [13, [l3 ' 
and entangled spins in electron gases due to exchange 
interaction [j, H, [13 ■ Usually entangling process takes 
two steps, generation of an entangled electron pair via 
some kind of interaction and separation of them from 
each other. 

A Cooper pair of a BCS superconductor is composed 
of two electrons with opposite momenta and a spin sin- 
glet state, k t and — k | [23|. Since the size of a Cooper 
pair is the order of the coherence length ^ [^ 10"'* cm), 
entanglement of electron spins may survive within that 
scale. This means the entanglement length may be about 
the coherence length. If one can extract a Cooper pair 
and separate two electrons from each other, the super- 
conductor may be a good natural resource of entangled 
spin states. In Refs. H H E El d ES El it is im- 
plicitly assumed that the distance between two tunnel 



junctions attached to the superconductor should be less 
than the coherence length ^. Also it is unclear whether 
a spin state of a tunneled electron pair is a spin singlet 
state (one of four Bell states) or a mixed state. 

In this Letter we address this problem and find the 
two-spin state of the BCS ground state is not given not 
by a Bell state but by a Werner state. We investigate 
entanglement of the two-spin state as a function of the 
relative distance between two electrons. Surprisingly, we 
find the entanglement length is not the order of the co- 
herence length ^ but the order of the Fermi wave length 
\f. 

Let us start with the pairing Hamiltonian of the BCS 
theory Ulllj 
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where c\^^ is a creation operator for electrons of wave 
vector k and z-component of spin s. The normalized 
BCS ground state of Eq. Q is given by 
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where coefficients Uk and Uk are written by 
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Here £'k — ^/(ek — /i)^ + A^ is the excitation energy of 
a quasi-particle of wave vector k and Ak the supercon- 
ducting gap. 

In order to investigate the entanglement of two- 
electron spins forming Cooper pairs, we introduce the 
two-electron space-spin density matrix in the second 
quantization 

p(2)(a;i,a;2;x;,4) = ]^{i>\x'^)i;\x\)^{x^)i>{x2)) , (4) 

where (• • • ) = (V'ol ■ ■ ' IV'o) at zero temperature and x — 
(r, s). The two-particle Green's function is defined by 

G(l, 2; 1', 2') = -(T[^H(l)^i/(2)V'l,(2')^lf (!')]> , (5) 



where 1 refers to Xiii, ■0}^ is the creation operator for 
electrons in the Heisenberg representation, and T the 
time-ordering operator. The relation between the two- 
electron space-spin density matrix and the two-particle 
Green's function reads 

p^'^\xi,x2;x'-i^,x'2) =--G{xiti,x2t2;x[t^,x'2t^), (6) 

where i^ denotes time infinitcsimally later than t. 

The two-particle Green's function for the supercon- 
ducting state can be factored into single-particle Green's 
functions ^E^ 



G(l, 2; 1', 2') = G(l, l')G(2, 2') - G(l, 2')G(2, l') 

-F(l,2)Ft(l',2'). (7) 

The single-particle Green's function is given by 

G(l, 1') = -z(T[7/^ff (xitO^Kt'i)]) (8a) 

^Ss,s',G{nh,r[t[), (8b) 

where its spin dependence for a non- magnetic system be- 



comes a unit matrix S^ 



The anomalous Green's func- 



tion is written by 

Ft(l,2) = ~t{T[^lf{^,h)i;Ux2t2)]) (9a) 

= Is,s,F\rih,r2t2) , (9b) 

where its spin dependence is given by an antisymmetric 
matrix 
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Also notice that F(l,2) = Is-^s2F{riti,r2t2). 

Since the system is translational-invariant and the 
Hamiltonian of the system is time-independent, the 
Green's functions G and F depend only on the relative 
coordinates ri — r[ and the time difference ii — 1[ . In the 
limit that t[ goes to ii, the spatial part of G becomes 



zG(r) = (V.o|V'i(ri)V^.(r;)|^o) 
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where 
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is the field operator in Schrodinger representation. In the 
continuum limit one obtains 



iG{r) 



1 m 

2^¥ JO 



Vf. &v[\(kr)k dk . 



(11) 



Due to the similarity between v^ of the BCS ground state 
and the Fermi function of an ideal electron gas at the 
critical temperature Tc, as shown in Fig. ^ Eq. I|ll|) has 



a similar form of Eq. (17) in Ref. M!|. The electron density 
n = N/V can be calculated by zG(0) = n/2. Similarly, 
the spatial part of F^ is given by 
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Note that F'l'(ri — r2) = [^(ri — ^2)]* ■ In the continuum 
limit one has 
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where ■fk = Ck — M, Ak — lS.9{fiWD — |^k|) with a step 
function 9{x), and the integration over k should be done 
on the range such that —Hlud < Ck < S^d |2(1 l2lll2a|. 
Here ojd is the Debye frequency. Thus Eq. (|13|l becomes 
approximately 

iFir)^Nm'^^M^) (14) 

where iV(0) is the density of states for one spin projection 
at the Fermi surface, given by [23| 
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and Ko{y) is a Bessel function of order 
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From Eqs. {7)), © and ©, we have the two-electron 
space-spin density matrix 

(2) / I l\ 

P.i,..;4,4(l"l'I"2;rir2) 

= -^[^..4'5..4G(ri-r;)G(r2-ri) 
-'5.,4<5,,,,G(ri-r^)G(r2-r;) 
- Is,sjs',,s', F(ri - r2)F*(r; - r^) ] . (17) 

Eq. p7(l has the same form of the two-electron space-spin 
density matrix for a non-interacting electron gas, Eq. (8) 
in Ref. [J], except the last anomalous term. In the limit 
that |ri — r^l — !■ 00 for i = 1, 2, one has 

(2) / I I \ 

Psi's2;ai,4(^l'r2;rir2) 

- \ls,sjs',.s', F{v, ~ r2)i^*(r; - r^ (18) 

which shows the off-diagonal long range order of a super- 
conductor |24| . 

Let us consider the case ri = r'j^, Y2 — r'j, which is 
equivalent to take only the diagonal elements of the space 
density matrix. For a solid state entangler that produce 
entangled spins out of a superconductor, two leads are 
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FIG. 1: (color online) v^ (solid line) and iikfk (dotted line) as 
a function of ek/eF. In plotting A/ep is taken to be 1/1000. 



attached to two tunneling points ri and r2 on the super- 
conductor. This implies one electron is located at ri and 
the other at r2. The two spin state depending on the 
relative distance r = ri — r2 is given by 
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Let us define functions g{r) = G(r)/G(0) — 2iG{r)/n 
and f{r) = F{r)/G{0). Eq. ^ becomes 
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where the two-spin state pi2 of the BCS ground state 
with normalization factor J\f = 4 — 2g^ + 2p is given by 
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We find that the two-spin state pi2 is not a spin sin- 
glet state but a Werner state characterized by a single 
parameter p {0 < p < 1) 

pi2 = (l-p)J+p|*(-))(*(-^|, (22) 

where I is a 4 x 4 unit matrix and |^(^^) = 
4= (I ti) — I it))- The parameter p is given by a function 
of / and g 
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For a non-interacting electron gas, the two spin-state is 
also given by a Werner state and can be obtained from 
Eq. inH) by putting / = 0. 

The properties of a Werner state are well known 4, 
|25| . According to the Peres-Horodecki separability cri- 
terion [23,123, a Werner state is entangled for p > 1/3. 
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FIG. 2: (color online) F (solid line), g (dotted line), C (thick 
solid line) as a function of kpr. The plotting parameters are 
given in the text. 



Thus P12 is entangled for p + 2g'^ > 1. If g^ = 0, pi2 is 
entangled for p > 1. As shown later, p is always much 
less than 1. The concurrence for pi2, one of entanglement 
measures, is calculated as 0, |23| 



C — max < 0, 



3p-l 



(24) 



Fig. 121 depicts the concurrence as a function of kpr. This 
implies that the two spin state pi2 of the BCS ground 
state is not entangled if the relative distance r between 
two tunneling points on the superconductor is larger than 
the Fermi wave length Xp and even though less than the 
coherence length f . 

Let us calculate the magnitude of f{r) which is rewrit- 
ten as 
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(25) 



where F{r) = F{r)/F{0) is always less then 1 except at 
the origin as shown in Fig. |21 On the other hand we 
obtain 
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In the weak coupling limit |20|, |21|, |23 , A < fuvo < ef, 
the value of Eq. (|26|l is very small. For example, for the 
values Awl meV, ?iajD w 100 meV, and ej? « 1 eV, we 
have F(0)/G(0) ~ 10^^. This means p is very small 
compared with g^. Thus we demonstrate that two spins 
are not entangled for Xp < r < £^. 

Our results could be explained as follows. A supercon- 
ductor has many Cooper pairs. Although the spin state 
of each Cooper pair is a Bell state, the spin correlation 
between two points ri and r2 is due to many Cooper pairs 
not a single Cooper pair. Thus for solid state entangler 
which generate entangled spins from a superconductor. 



one needs to extract a single Cooper pair or to confirm 
that two electrons are tunneled out from a single Cooper 
pair. 

In conclusion, we obtained the two-electron reduced 
density matrix of the BCS ground state based on the 
Green's function method. We investigated entanglement 
of two electron-spins forming the Cooper pair in BCS 
superconductors. It has been found that the two-spin 
density matrix for a given relative distance between 
two electrons is given by a Werner state not by a Bell 



state. Also the entanglement length is not the order of 
the coherence length ^ but the order of the Fermi wave 
length A_F. 
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